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It has by now been established that standard QCD factorization using transverse momentum de- 
pendent parton distribution functions fails in hadro-production of nearly back-to-back hadrons with 
high transverse momentum. The essential problem is that gauge invariant transverse momentum 
dependent parton distribution functions cannot be defined with process-independent Wilson line op- 
erators, thus implying a breakdown of universality. This has led naturally to proposals that a correct 
approach is to instead use a type of "generalized" transverse momentum dependent factorization 
in which the basic factorized structure is assumed to remain valid, but with transverse momentum 
dependent parton distribution functions that contain non-standard, process dependent Wilson line 
structures. In other words, to recover a factorization formula, it has become common to assume 
that it is sufficient to simply modify the Wilson lines in the parton correlation functions for each 
separate hadron. In this paper, we will illustrate by direct counter-example that this is not possible 
in a non-Abelian gauge theory. Since a proof of generalized transverse momentum dependent fac- 
torization should apply generally to any hard hadro-production process, a single counter-example 
suffices to show that a general proof does not exist. Therefore, to make the counter-argument clear 
and explicit, we illustrate with a specific calculation for a double spin asymmetry in a spectator 
model with a non-Abelian gauge field. The observed breakdown of generalized transverse momen- 
tum dependent factorization challenges the notion that the role of parton transverse momentum in 
such processes can be described using separate correlation functions for each external hadron. 
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I. INTRODUCTION 

There has been much recent activity devoted to the 
study of parton transverse momentum in high energy 
hadronic collisions. Observables that are sensitive to par- 
ton transverse momentum can potentially provide new 
insight into the structure of hadrons. However, to in- 
terpret data it is important to understand the extent to 
which familiar perturbative QCD (pQCD) factorization 
approaches can be extended to situations that involve 
parton transverse momentum. In this paper, we specif- 
ically address recent efforts to apply the usual general 
framework of pQCD factorization (with some modifica- 
tions) to describe parton transverse momentum in colli- 
sions between high energy hadrons with production of a 
pair of nearly back-to-back high transverse momentum 
hadrons or jets in the final state: 

H l +H 2 ^H 3 + H i + X. (1) 

The distribution of transverse momentum inside the col- 
liding hadrons can be probed, for example, by measuring 
the small imbalance in the distribution of transverse mo- 
mentum between the final state pair (see, e.g., [H, HI)- 
Parton transverse momentum also plays a central role in 
the generation of spin asymmetries. 
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For a discussion of the relevant issues, it will be impor- 
tant to first clearly define our terminology. Very gener- 
ally, a QCD factorization theorem [3[ is said to be valid 
if an observable (such as a cross section) can be written 
as a convolution product of factors that describe different 
regions of parton momentum. Schematically, one expects 
for process ([T]): 

dcr = n (g) (g) $ ff2 (g) A Hs <S> Aff 4 + p.s.c. (2) 

Here, the symbol (g> denotes all relevant convolution in- 
tegrals and traces, and a sum over different types of par- 
tons and subprocesses is understood. The hard factor H 
describes the short range behavior in the hard collision 
between partons, while and <&h 2 are the parton dis- 
tribution functions (PDFs) for initial state hadrons Hi 
and i?2- The A# 3 and A# 4 are fragmentation functions 
for final state hadrons H3 and H4. In general, a soft 
factor may also be needed, though we will not write it 
explicitly. The symbol "p.s.c." indicates that power sup- 
pressed corrections are neglected. 

A transverse momentum dependent (TMD) factoriza- 
tion theorem (also called unintegrated factorization or 
fcT-factorization) is said to be valid if the role of parton 
transverse momentum can be taken into account by us- 
ing PDFs and FFs in Eq. © which depend explicitly on 
parton transverse momentum. In a TMD-factorization 
formula, one therefore refers to TMD PDFs and TMD 
FFs (or "unintegrated" PDFs and FFs). 

TMD-factorization should be contrasted with the 
more common collinear factorization theorems, applica- 
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ble to cases where observables are not sensitive to in- 
trinsic transverse parton momentum. In the collinear 
factorization theorems, transverse momentum is inte- 
grated over inside the definitions of the standard PDFs 
and FFs. These are the standard "integrated" PDFs 
and FFs which are expressible as well-defined, process- 
independent operator matrix elements Q. The process- 
independence means that the integrated PDFs and FFs 
can be parameterized by experimental data and then 
later reused in calculations to make first-principle pre- 
dictions for future experiments. Thus, this universal- 
ity property of the standard integrated PDFs and FFs 
lends great predictive power to the standard pQCD ap- 
proaches, and is a basic component of the standard 
collinear factorization theorems. 

It is natural to hope that an analogous universality 
property applies to the TMD PDFs and TMD FFs in 
a TMD-factorization formula. In that case, we would 
say that standard TMD-factorization is valid, where we 
will explicitly use the word "standard" to refer to the 
universality condition. If a standard TMD-factorization 
formula were valid for the process in ([1]), then it would 
be possible to calculate cross sections using the same 
TMD PDFs and FFs that are parameterized in other 
processes like deep inelastic scattering (DIS) or the Drell- 
Yan (DY) process. However, as we will discuss in more 
detail, TMD-factorization is not generally valid for the 
process in Eq. ([1]). 

The essential complication arises from longitudinally 
polarized gluons that couple soft and collinear subgraphs 
to the hard part and which, at first sight, appear to 
break topological factorization graph-by-graph at lead- 
ing power. Dealing with these "extra" gluons is one of 
the main issues that must be dealt with in all factoriza- 
tion proofs in pQCD. In DIS, for example, it is found that 
after summing over all graphs (and applying appropriate 
approximations), a Ward identity argument allows the 
extra gluon contributions to be factored into contribu- 
tions which correspond to Wilson lines (also called gauge 
links) in the definitions of the PDFs or FFs 3]. How- 
ever, in order to justify the approximations that allow an 
application of the Ward identity, it is necessary to first 
perform certain contour deformations on the gluon mo- 
mentum For DIS, the contour deformations must be 
consistent with having extra gluon attachments between 
hadron spectators and a final state struck quark. Similar 
steps apply to DY, but there the contour deformations 
should correspond to extra attachments between hadron 
spectators and an initial state quark. (See also Ref. @.) 
The result is that the Wilson line for the quark TMD 
PDF in DIS is future pointing while the Wilson line for 
a quark TMD PDF in DY is past pointing Q. The dif- 
ference in Wilson line direction has been shown to result 
in a sign flip for the Sivers function (a particular type 
of TMD PDF) in DY as compared to DIS 0. Strictly 
speaking, this could be regarded as a breakdown of stan- 
dard TMD-factorization because the Sivers function is 
not truly universal in DIS and DY. However, since a sign 



flip is easily accounted for, it is more appropriate to say 
that the Sivers function in DIS and DY possesses a type 
of "modified" universality. The relationship between cor- 
relation functions in DIS and DY is further discussed in 
Ref. 3. 

The problems that occur with TMD-factorization for 
Eq. ([1} are much more complicated because there one 
must deal with extra gluons that connect spectators 
to both initial and final state partons. The result is 
that the contour deformations necessary for a factor- 
ization proof prevent a direct application of the usual 
Ward identity arguments. This problem was observed 
by Bomhof, Mulders and Pijlman [9j who found that the 
Wilson lines needed for gauge-invariant TMD PDFs and 
FFs in hadro-production of hadrons are not generally 
process- independent. Although the extra gluon attach- 
ments eikonalize, the resulting sums of eikonal factors 
do not correspond to the simple future or past pointing 
Wilson lines that are found in DIS and DY. In [9( it was 
shown at the level of a single extra gluon that the TMD 
PDFs and FFs (assuming consistent definitions exist) are 
non-universal because they require, at a minimum, pro- 
cess dependent Wilson lines. Hence, there is a violation 
of the universality property necessary for standard TMD- 
factorization. 

Similar problems are encountered in proofs of collinear 
factorization for hadron-hadron scattering, but there one 
is saved by cancellations between graphs that occur after 
integration over transverse momentum (see, e.g., Ref. Q 
and also more recent work in Ref. [10j ) . Since these can- 
cellations are not point-by-point in transverse momen- 
tum, they do not generally apply to TMD-factorization. 

To address the role of non-universal Wilson lines, the 
concept of a "generalized" TMD-factorization formula 
was later developed |ll|, [l2|. In this approach, it is 
assumed that the only deviation from standard TMD- 
factorization is that the TMD parton correlation func- 
tions (PDFs and FFs) must contain process-dependent 
Wilson lines. Schematically, instead of Eq. © one must 
assume a more complicated expression of the form 

da = 

+ p.s.c. (3) 

Here, H c ,j is the hard part for subprocess j and color 
routing c. For each subprocess and routing of color 
through the hard part, there is in general a different set 
of TMD PDFs and TMD FFs corresponding to the dif- 
ferent Wilson line structures that they contain. In the 

above notation ( i > ^^ 1 ' is a gauge invariant TMD PDF 

for hadron Hi with a Wilson line [WLl' 3 ] correspond- 
ing to subprocess j and color routing c. Analogous no- 
tation is used for the other process-dependent correla- 
tion functions. (There are also possible soft factors not 
shown explicitly in Eq. Each term in Eq. © has 
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the same basic factorized structure as in Eq. ((2J, involv- 
ing distinct (though process-dependent) TMD PDFs and 
FFs for all external hadrons. The TMD PDFs and FFs 
have the usual structure of a pair of field operators and 
a Wilson line with an expectation value corresponding 
to a specific external hadron state. The only difference 
from the standard case is that they are equipped with 
non-standard and potentially complex Wilson line struc- 
tures. In particular, Eq. (|3]) contains no matrix element 
of the form {H l H 2 \ ■ ■ ■ \Hi_H 2 ). So, by generalized TMD- 
factorization we mean that a TMD-factorization formula 
is recovered simply by replacing the Wilson lines in the 
definitions of the correlation functions by non-standard 
ones, which may be different for each hard subprocess and 
for each way of routing color through the hard part. The 
TMD PDFs needed in a generalized TMD-factorization 
formula for Eq. (UJ) could be totally different from the 
ones parameterized in, e.g. DIS and DY. 

A conjectured TMD-factorization of the form of Eq. J3| 
is a basic assumption in a number of recent studies [lll - 
IT^ |. The minimal Wilson line structures needed for 
Eq. ([3]) can be determined by considering a single ex- 
tra gluon at a time, radiated from each of the ex- 
ternal hadrons and attaching everywhere in the hard 
subprocess. The resulting process-dependent gauge in- 
variant correlation functions have been tabulated in 
Refs. [HI, [HI- These correlation functions have also been 
used to calculate physical observables such as weighted 
spin asymmetries [T3I . 

Collins and Qiu [ItJ verified explicitly that standard 
TMD-factorization fails in a sample calculation of a sin- 
gle spin asymmetry (SSA). That is, they showed in an ex- 
plicit calculation that the process-dependence of the Wil- 
son line structures observed in Ref. Q indeed corresponds 
to non-universality for the TMD PDFs. For their calcu- 
lation, they used a model Abelian theory and calculated 
the effect of a single extra gluon. An explicit illustration 
of the violation of standard TMD-factorization was also 
given for unpolarized scattering in Ref. |la |. again using 
a model Abelian gauge theory. The two-gluon example 
for unpolarized scattering is important as it directly il- 
lustrates that standard TMD-factorization cannot gener- 
ally be recovered by rescaling the hard part with a con- 
stant color factor. (Co mpare this with the procedure of 
Refs. [l9|.) In Ref. [2Q|. it was shown explicitly that the 
observed breakdown of standard TMD-factorization de- 
scribed in Refs. [13, [3 is consistent with the gene ralized 
TMD-factorization proposed in Refs. [Til [T^. ITEj. again 
within the Abelian theory. 

However, all the cases studied so far have only con- 
sidered graphs with extra gluons radiated from one of 
the hadrons at a time. What is missing is a treatment 
of non-Abelian gluons radiated from different hadrons 
simultaneously. If a generalized TMD-factorization ap- 
proach is possible, then extra gluons radiated from all 
hadrons simultaneously must be shown to eikonalize and 
factorize after a sum over graphs. Given the complex 
color structures that arise in a non-Abelian gauge the- 



ory, it is unclear that such a procedure is possible in real 
QCD. 

The purpose of this paper is to show explicitly that 
even generalized TMD-factorization breaks down in a 
non-Abelian gauge theory at the level of two extra glu- 
ons. In other words, the violation of standard TMD- 
factorization, already found in previous work, cannot be 
dealt with simply by replacing the Wilson lines in the 
standard correlation functions by more complicated ones 
and summing over different subprocesses and color struc- 
tures as in Eq. ©. 

As seen in Ref. 
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down of standard TMD-factorization are illustrated most 
directly in a calculation of an SSA with a single extra 
gluon. We will find analogously that the breakdown of 
generalized TMD-factorization in a non-Abelian gauge 
theory is most easily illustrated in a calculation of a dou- 
ble Sivers effect in a double transverse spin asymmetry 
(DSA). As in Refs. [HEl we will use a model field the- 
ory to describe the quarks, spectators and hadrons. A 
proper counterexample to generalized TMD-factorization 
must verify that terms which violate generalized TMD- 
factorization graph-by-graph do not cancel in a sum over 
graphs. This is most easily done in a simple spectator 
model that restricts the number of relevant Feynman 
graphs. 

In Sect. [TT] we discuss the particular model and de- 
scribe the procedure for deriving a violation of gener- 
alized TMD-factorization. In Sect. QTlJ we review the 
steps for factorization with one extra gluon. We explic- 
itly review the breakdown of standard TMD-factorization 
for two extra gluons from one hadron in Sect. IIV1 In 
Sect.[V]we discuss the generalized TMD-factorization for- 
mula that is required to recover a factorized structure. 
In Sect. IVII we demonstrate that the generalized TMD- 
factorization formula is inconsistent with having extra 
gluons radiated from both hadrons simultaneously. We 
end with concluding remarks in Sect. IVTll 



II. SETUP 

A simple model field theory provides a direct illustra- 
tion of why factorization fails in a gauge theory, while 
avoiding the complications of dealing with a large num- 
ber of Feynman graphs. We will continue to use the 
model field theory of Refs. 0, H H, though with a 
few important differences. The hadrons continue to cor- 
respond to different flavors. The "quarks" continue to 
correspond to scalar fields </>/, while the "hadron" fields 
Hf and the spectator "diquarks" ipf are Dirac spinors. 
The subscri pt f = 1, 2 labels flavor. The main difference 
from Refs. HI [H, [2l| which we will introduce is that 
the gauge field will be the massless SU(7V C ) non-Abelian 
gauge field. (In QCD N c = 3.) By contrast, Refs. [13, [3 
used a massive Abelian gauge field that coupled with dif- 
ferent charges, g\ and g%, to the quarks and diquarks 
in hadrons H\ and H2- In this paper, the non-Abelian 
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FIG. 1: Basic graph contributing to hadro-production with no extra gluons. 



gauge field couples with equal strength to the quarks in 
each hadron, g\= gi = g. 

We also will introduce one more field: The hard sub- 
process will be described by an exchange of a hypotheti- 
cal new massive color-neutral U(l) gauge boson which we 
will call X . It couples with strengths A^ 2 ) to quarks of 
flavor 1(2). Hence, the hard subgraph has a trivial color 
structure. This will allow us to study the role of extra 
gluons in a non-Abelian gauge theory, while limiting the 
number of Feynman graphs that need to be explicitly 
considered, and eliminating color flow directly through 
the hard part. It is important to emphasize that the 
general conclusions of this paper do not depend on this 
specific model for the hard subprocess. What is impor- 
tant is that color is carried by both the initial and final 
state partons. A proof of generalized TMD-factorization 
for such processes should be general and apply to any 
hard hadro-production process with observed hadrons in 
both the initial and final states. If it fails for a specific ex- 
ample, then no general proof exists. In cases where color 
is exchanged in the hard process, the Wilson line struc- 
tures for the TMD PDFs in Eq. © are determined after 
first expanding the color matrices of the hard subprocess 
into different routings following the procedure described 
in Ref. [121 ] . The steps in this paper then apply to each 
term in a series of color routings, and similar contradic- 



tions with generalized TMD-factorization can be found. 
Our choice of hard subprocess is to make the illustra- 
tion of why generalized TMD-factorization fails as clear 
as possible. 

The basic graph contributing to hadro-production of 
high-p t hadrons, with no extra gluon attachments, is 
shown in Fig. [TJ Because of the trivial color structure of 
the hard part, and the simplicity of our model of quarks 
and spectators, there is only one subprocess and one color 
flow. So a generalized TMD-factorization formula corre- 
sponding to Eq. ([3|) can only involve one term. 

We work in the center-of-mass frame where the incom- 
ing hadron H\ is initially moving with large rapidity in 
the forward plus direction with no transverse momen- 
tum, while H 2 is initially moving with large rapidity in 
the minus direction with no transverse momentum. Each 
hadron splits into an active "quark" which enters the 
hard subgraph, and a spectator "diquark" which enters 
the final state. The quarks interact in the hard part by 
exchanging a hard colorless vector boson X, with momen- 
tum q where \q 2 \ is large, |g 2 | ^> M\. Within the model, 
the final state high transverse momentum hadrons or jets 
are represented simply by on-shell final state quarks. 

The polarization dependent differential cross section at 
zeroth order in g is, 
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(tf -m 2 qi ) 2 



,2 \2 



(^2 ""<J2^ 

x (2ir) 2 6(( Pl - fci) 2 - m^)«((p2 



k 2 ? 



m% 2 ). (4) 



The factor of N 2 comes from tracing over the unit color triplet matrices in the upper and lower color loops in Fig. [T] 
We will use the approximation that the struck partons are collinear to the directions of their parent hadrons. So, 
k^ ~ and k 2 ~ p 2 . The minus component of k\ and the plus component of k 2 are small (order A 2 /p^) and can be 
neglected inside the hard part. We also define x\ = k\ jp\ and x 2 = k 2 /p 2 , and s is the usual Mandelstam variable 
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FIG. 2: Typical cases of graphs that vanish when extra gluons are considered because of the trivial color factor, Trc [t a ] = 0. 



for center-of-mass energy squared. The transverse spin vectors Si( 2 ) label the transverse spin for hadrons #1(2) an d 
are normalized such that the extreme values are = 1. The subscript D on Ttd [• ■ •] indicates a trace over Dirac 
indices. 

After evaluating the (^-functions, this cross section easily factorizes: 
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d 3 k 3 d 3 k 4 2s J (2tt) 2 



16tt 3 



k 2 

"'IT 



xi(l - xi)m 2 Hi + (1 - a;i)m^ + scim^ 



7V c A2x 2 (l - x 2 ) 2 



^Tr D 



+ mH 2 )(l + 75? 



^ 2 + "V 2 ) 



where the hard part is, 



16tt 3 



A; 2 A 2 2 f (fci + fc 3 ) • (fcs + 2fc 4 - fci) 



x 2 (l - x 2 )m 2 H2 + (1 - x 2 )m2 2 + X2 



H(ki,k 3 , fc 4 ) 



2xix 2 s 



(fcj - fc 3 ) 2 - Af 



x 



(5) 



(6) 



The remaining factors in braces in Eq. ([5]) are what is expected from the zeroth order expansion of the standard 
operator definition of a TMD PDF. For hadron Hi, for example, the zeroth order TMD PDF is 



0(1) _ iVcAfgi(l ~ x i) 



±Tr T 



16tt 3 



1,2 
"'IT 



• £l(l - Xi)m 2 Hi + (1 - x^m 2 ^ + x\i 



(7) 



and similarly for hadron H 2 . Note that k 2 T is not kine- 
matically independent of \tyr- k 2 T = k 3 T + k 4 T — kxy. 
Also, x\ » (fcg + and £ 2 w (fcg + k^)/p 2 . To 

keep the expression simple, we have written Eq. ([5]) in 
terms of exact parton momenta, fci,fc 2 ,fc 3 ,fc 4 , though we 
remark that in the correct final formula these should be 
replaced by the approximate on-shell parton momenta 
appropriate for a hard subgraph |3l| . 

When treating higher order graphs, we will always 
work in Feynman gauge, where the contour deformations 



needed in a derivation of factorization are most straight- 
forward. For factorization to work, it must be possible to 
identify any uncanceled collinear or soft singularities as 
contributions to non-perturbative correlation functions. 
Some of these singularities correspond to Wilson line con- 
tributions, and are therefore essential for maintaining 
gauge invariance. 

Before continuing we should mention that there are 
a number of general complications involved in deriving 
TMD-factorization that will not be addressed here be- 
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cause they are not directly related to the main reasons 
that generalized TMD-factorization fails. The most naive 
definitions of TMD correlation functions include extra 
divergences that need to be removed by appropriate re- 
definitions (see Ref. [22| and references therein for an 
overview of issues related to the precise definition of a 
TMD PDF). One particular complication is that the Wil- 
son lines in these definitions cannot be exactly light-like 
without containing extra "light-cone" divergences which 
correspond to partons moving with large rapidity oppo- 
site to the direction of their parent hadrons. Proposed 
solutions involve either tilting the Wilson line away from 
the exactly light-like direction USL or dividing out by ex- 
tra gauge invariant factors 0, l23l - [25| . Furthermore, an 
exactly correct factorization formula requires a soft fac- 
tor to account for gluons with all components small. The 
other parton correlation functions tend to overlap in the 
soft region, and a fully correct definition of a TMD PDF 
requires extra factors to remove the overlap. 

There has been much significant recent work devoted 
to finding a fully consistent definition of the TMD PDF. 
However, our main concern in this paper is only with 
the general color structure of the Wilson lines inside the 
matrix elements for the external hadrons. Therefore, for 
our purposes it will be sufficient to continue to treat the 
direction of the main Wilson lines in the TMD PDFs as 
being light-like, and we will not address the role of a soft 
factor or the overlap of regions. Finally, we will restrict 
consideration to the limit of very large relative transverse 
momentum (large-pt) where one expects standard pQCD 
methods to be most appropriate, and we do not consider 
the possibility of recovering a type of TMD-factorization 
appropriate in the small- a; limit 26] . 

Having discussed the basic graph in Fig. [TJ the next 
step is to consider graphs dressed with extra gluons. The 
graphs which can contribute to the Wilson line insertions 
are those in which extra, nearly on-shell gluons connect 
different subgraphs — for example, graphs with gluons 
connecting the pi-collinear lines to the outgoing struck 
quark lines or to the p2-collinear lines. Normally, one 
expects the sum of such graphs to contribute to the Wil- 
son line in the TMD PDF for Hi, after application of 
a Ward identity. The primary issue is that the approx- 
imations that normally allow Ward identities to be ap- 
plied are only valid after certain contour deformations 
on the extra momentum integrals. Namely, to apply a 
Ward identity, it must be possible to approximate an ex- 
tra gluon by a longitudinally polarized one with a large 
component of longitudinal momentum. Then the extra 
gluon momentum can be contracted with the hard scat- 
tering matrix element and a Ward identity argument can 
be applied directly. However, when virtual gluons at- 
tach to a spectator line, they give contributions from 
the "Glauber" region, meaning that if / is the momen- 
tum of an extra gluon, then \l + \ and |/ _ | are both much 
smaller than \It\- In the Glauber region, the approxi- 
mations needed for the Ward identity are not valid. For 
factorization to work, it must be possible to first deform 



the contour out of the Glauber region (see, for example, 
Refs. [(|). Alternatively, factorization could be recovered 
if there is a cancellation between graphs, as in the stan- 
dard proofs of integrated (collinear) factorization. The 
basic pro blem with TMD-factorization found in previous 
work [TP . [IH [TH, [13, EH is that the necessary contour 
deformations needed to treat the Glauber region are in- 
consistent with a direct application of the standard Ward 
identity arguments. Namely, they are in different direc- 
tions for different graphs depending on whether the in- 
teraction is in the initial or final state. 

So, we will only consider graphs that can yield con- 
tributions from the Glauber region. Since real gluons 
can never be in the Glauber region, we will only con- 
sider graphs with virtual gluons. Also, as long as no 
restrictions are placed on the target remnant momenta, 
graphs with spectator-spectator interactions cancel |5| in 
the integration over final states [32[. Similar cancella- 
tions occur between different cuts of the same graph for 
active-spectator interactions after parton transverse mo- 
mentum is integrated over, and are needed in the stan- 
dard proofs of collinear factorization [j| . A counter- proof 
of TMD-factorizaton therefore needs to show that such 
cancellations generally fail when transverse parton mo- 
mentum is explicitly taken into account. A specific ex- 
ample of such a non-cancellation was given in Ref. (l8j 
and will be reviewed in Sect. IIVI In graphs with attach- 
ments between active quarks, there are not enough Dirac 
7-matrices to give spin dependence to the TMD PDFs. 
Such graphs will therefore not affect our discussion of sin- 
gle and double spin asymmetries at lowest non- vanishing 
order. Furthermore, graphs with a scalar-scalar-gluon- 
gluon vertex do not give leading power contributions to 
cikonal factors. 

We remark that, because the TMD factorization break- 
ing effects are due to the Glauber region where all com- 
ponents of gluon momentum are small, the interactions 
responsible for breaking TMD factorization are associ- 
ated with large distance scales. 

In our specific model, a large number of graphs van- 
ish simply because of the highly simple color structure 
involved. Examples are shown in Fig. [5J They vanish 
because their color factors include a trace around a color 
loop of a single SU(iV c ) generator, Tr c [t a ] = 0. (The C 
on the Trc [■ ■ '] denotes a trace over triplet color indices.) 

Hence, the relevant types of graphs are represented by 
Figs. [3] through [8] If a generalized TMD-factorization 
formula is possible, then the sum over all such graphs 
must produce a factorized form like Eq. §S§ with a Wilson 
line structure in the TMD PDF or FF for each hadron 
separately. We will consider each type of graph in the 
following sections. 



III. ONE EXTRA GLUON 

We begin the investigation of diagrams by reviewing 
the steps for determining the contribution from a single 
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FIG. 3: Extra gluon attachments that are consistent with the standard Wilson line structure. The difference in shading on 
the arrow at the eikonal attachments for (a,b) and (c,d) is to emphasize that the gluons are from different hadrons. They 
correspond to pi in (a,b) and to P2 in (c,d). 



extra gluon. As in Ref. [17]], we focus on the calculation 
of an SSA. We start with graphs of the type shown in 
Fig. [3j where the extra gluon attaches on the side of the 
hard part nearest to its parent hadron. Any spin asym- 
metry disappears in the zeroth order cross section, Fig. [1] 
because there are too few Dirac matrices to produce a 
non-zero result in the traces with 75. 

Consider, for example, Fig. [21a). The arrow on the 
gluon line indicates that it is collinear to Hi. By first 
deforming the I integral out of the Glauber region to 
the iii-collinear region, one may replace the intermediate 
struck quark line of momentum k^—lhy the eikonal factor 



t a gn^ 
-l+ + ie 



-gt a n^PN.— - igt a n^TrS(l+) (8) 

where = (0, 1, t ). The sign on the ie is determined by 
the direction of the contour deformation. For the spin- 
dependent part, the attachment of the extra gluon at the 
spectator produces a factor at leading power equal to 



t" , 



-Tr D 



i> 1 + niH l )l5^ 1 {f 1 - %i + t + WtyJ x 
w 2it a e 3k s[l k p + (m Hl (l - xi) + m 01 ). (9) 



When this expression is combined with the imaginary 
part of Eq. ©, the factors of —i and i combine and a 
contribution to an SSA is obtained. The €jk is the two- 
dimensional Levi-Civita symbol with ei2 = 1. 

If the extra gluon is on the other side of the cut as in 
Fig.GHb), the eikonal factor is 

= -gt a n>£P.V± + ig t a n^S(l+). (10) 
The factor from the attachment at the spectator is, 



t a 

-Tr D 
2 D 



\ + ™hJ75^ 1 (|> 1 - h + m^ x ) x 

x 7 + (| > 1 -#!+/+ "tyi 

« -2it a e jk s{l k p + {m Hl {l -xi) +ra^). (11) 

When this is combined with the imaginary part of 
Eq. (fTU)) , the factors of i and —i combine and a contribu- 
tion to an SSA is again obtained. It is exactly equal to 
the one found in Fig. [3ta) as it must be since the graphs 
are related by Hermitian conjugation. 

Analogous steps apply to Figs. [2Ic,d). The different 
shadings on the arrows is to emphasize that they are ra- 
diated from H2 rather than H\. Now, after a contour 
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FIG. 4: Graphs of the type that led to a violation of TMD-factorization in Ref. [lj- Gluons radiated from Hi are illustrated 
in (a) and (b); gluons radiated from H2 are illustrated in (c) and (d). The Hermitian conjugate graphs should also be included. 
These graphs vanish in our model because of their trivial color factor. 



deformation I can be made collinear to the minus direc- 
tion, and a vector n 2 = 
The eikonal factors are, 



(1,0, Ot) is used instead of 



igt a n 2 ir 5(1 



(12) 



-gt a n%P.V.— + igt a n v 2 ir5{l-) (13) 



-l~ + ie 
for Fig. rj^c) and 

t a gn» 2 = 
—l~ — ie 
for Fig.^d). 

In our calculation of the SSA, p 2 is unpolarized so there 
is no factor of i coming from the attachment to the up- 
per spectator. Therefore, it is only the real principal 
value contributions that are kept in Eqs. (|12I13[) . Note 
that there can be no double Sivers effect at the level of 
one extra gluon because there must be at least one extra 
gluon from each of the hadrons for their PDFs to have 
spin dependence. 

The t a s in Eqs. (|8ll0p combine with the t a s in 
Eq. (|9llip to produce an overall color factor equal to 



T F {N 2 C - 1). 



(14) 



Exactly the same color factor also arises when I is 
collinear to p 2 as in Figs. GHc,d). 



The remaining type of single-gluon graphs are those in 
Fig. |U For the lines collinear to H\ the eikonal propaga- 
tors are 

Tr c [t a ] gr4 = 
-1+ + ie 

- 5 Tr c [t a ] <P.V.^ - tgTvc [t a ] n^S(l+), (15) 

from Fig. |4ja) and 

Tr c [t a ] gr4 _ 
1+ + ie 

+ gTr c [t a ]n^P.V.^-igTr c [t a ]n^S(l + ), (16) 

from Fig. HJb). The sum of these graphs gives a contri- 
bution equal to 



- 2igTr c [t a ] 7T<<5(Z+) 



(17) 



If the color factor in Eq. (fTTl) did not vanish, then it would 
contribute to an SSA when the factor of i combines with 
the factor of i in Eq. ©. The same result is obtained 
from the Hermitian conjugate graphs. In a totally un- 
polarized cross section, there would be no contribution 



9 



from Eq. (IT7|) . regardless of the color factor, since it is 
imaginary. 

A symmetric analysis applies if the extra gluon is in- 
stead radiated from hadron H 2 as in Figs.2tc,d). In that 
case rii is replaced by n 2 and l + is replaced by l~. Since 
p 2 is unpolarizcd, it is the real parts of the eikonal fac- 
tors that contribute. Thus, since the real parts cancel 
between graphs, the eikonal factors from Fig. 0Jc,d) do 
not contribute to the unpolarized TMD PDF of hadron 

The uncanceled terms like Eq. p7|) would ordinarily 
signal a breakdown of standard TMD-factorization be- 
cause they are not consistent with the standard Wil- 
son line structure in a TMD PDF. This is exactly what 
is observed for the Abelian gauge theory calculation 
in Ref. [17j . In our non- Abelian example, however, 
these single extra gluon contributions to a "factorization 
anomaly" are exactly zero because they happen to in- 
clude the trivial color factor: 



Tr c [t a ] = 0. 



(18) 



So, in our specific non- Abelian calculation, Eqs. (0 [TU)) 
are the only eikonal propagators that contribute to an 
SSA at the level of just one extra gluon. 

The non-vanishing eikonal factors (coming from 
Figs. |3Ja,b)) are exactly what is obtained from an order 
g expansion of the Wilson lines in the standard definition 
of the TMD PDFs: 



dw <i 2 w t 



-iiiPj tu +ik t -w t . 



x {Hx , Sl \4 d (0, uT, w t ) U% l] [0, w] M (0) Iff! , 8l ) . 



(19) 



The standard Wilson line operator [0, 10] is inserted 
between the two scalar quark fields. The triplet color 
indices j, k are shown explicitly to emphasize the flow of 
color. The full Wilson line insertion is, 

u$ l] [o,w] = [Viim)].., [i{ ni )] rk , [V (m)} k , k , (20) 

where 

[V w (ni)] n , = 

Pexp(-igt a d\ ni ■ A a (w + A«i) J , (21) 
V Jo Jjf 

with P denoting a path-ordering operator. The super- 
script [rii] in Eq. (|19[) refers to the direction of the main 
leg of the Wilson line, starting from point 0. The ex- 
tra gluon attachment in Fig. [3^a) on the left side of the 
cut contributes to [Vo(ni)] k , k while the extra gluon in 
Fig. fflb) contributes to [V^(ni)].. r . To close the Wil- 
son line and ensure that it is exactly gauge invariant, the 
path needs a transverse link at light-cone infinity [271 ] . 



[I(ni)y = Fexp (-igt a J dz'Afr 



where C is a path in the transverse direction connect- 
ing the points, (0,oo,w t ) and (0,oo,O t ). In a deriva- 
tion of TMD-factorization in Feynman gauge, contribu- 
tions to the link at infinity do not arise explicitly. (See, 
e.g., Refs. SH and references therein for a review of 
the steps for resumming collincar gluon attachments and 
identifying the resulting gauge links.) 

Note that the Wilson line insertion in the standard 
TMD PDF of Eq. (JTSJ contracts the color indices of the 

quark fields, so {x\ , kvr) has no leftover color indices 
— it is a real- valued function as is appropriate for a gauge 
invariant probability density. 

Similar steps result in the single-gluon contribution to 
the Wilson line in the standard TMD PDF for H 2 when 
the extra gluon is radiated from the upper spectator as 
in Figs. G2c,d). There the Wilson line insertion instead 
points in the direction n 2 = (1,0, Ot): 



dw + d 2 vs t ZX2P - w + +zkt . Wt . 
(2tt) 3 



(22) 



3 j 



x (H a , s 2 \4 . (0, w + , w t ) U\ n k 2] [0, w] ^ k (0) \H 2 ,s 2 ). 

(23) 

In the Abelian case considered in Ref. breaking 
of standard TMD-factorization was observed for an SSA 
because the factorization anomaly terms analogous to 
Eq. (|17l) did not vanish. There, the "color factor" was the 
Abelian charge g 2 for Abelian "gluons" rather than the 
non-Abelian color factor Trc [t a ] = 0. That the anoma- 
lous terms vanish at the level of one extra gluon in our 
non-Abelian example is only due to the highly simple 
color structure in the particular hard process that wc 
have considered. Generally, when color is exchanged in 
the hard part, standard TMD-factorization breaking will 
already appear at the level of one extra gluon. In our 
example, violations of standard TMD-factorization only 
appear at the level of two extra gluons or higher. 

Schematically, the standard TMD-factorization for- 
mula suggested by the sum of one-extra-gluon graphs in 
Figs, ifa-d) is 

da = H® ^ (m, fci T ) ® $£ 2 21 fa, k 2T )® 

® <S (2) (ki T + k 2T - k 3T - k 4T ) . (24) 

Here H is the same hard factor that appeared at zeroth 
order in Sect. [IT] The "!" on the equal sign is to emphasize 
that this formula is ultimately incorrect, as explained in 
Refs. pj], [3 • We will see this explicitly at the level two 
extra gluons in the next section. 



IV. TWO EXTRA GLUONS 

In this section we review the steps for illustrating a 
violation of standard TMD-factorization at the level of 
two extra gluons radiated from one of the hadrons. The 
structure of the resulting TMD-factorization anomaly 
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FIG. 5: A gluon three-point interaction that contributes to the single gluon Wilson line attachment. 
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FIG. 6: Graphs that contribute to the standard Wilson lines at the level of two extra gluons. The Hermitian conjugate graphs 
should also be included, as well as graphs with both gluons collinear to p2- 



terms constrains the possible generalizations of TMD- 
factorization, as we will discuss more in Sect. [V] 

The simple color structure of the hard process strongly 
limits the number of non- vanishing graphs. All graphs in 
which the second gluon is internal to the hadron sub- 
graph (e.g. Fig. [S]) only give further contributions to a 
single gluon contribution in the standard Wilson lines in 
Eq. (|24[) . The only non- vanishing graphs that give order 
g 2 contributions to a Wilson line operator are of the type 
shown in Figs. [51 El [HI All other graphs vanish because 



they include an overall color factor Trc [t a ] = 0. 

We will consider each type of graph in turn in the 
next few subsections. We will first identify the graphs 
that are consistent with the standard TMD-factorization 
formula Eq. (|24l) . Next we will consider the graphs that 
contribute to a violation of standard TMD-factorization 
that were already discussed in Refs. (l8| . The remaining 
graphs in Fig. [5] are ultimately responsible for breaking 
generalized factorization, so we will put off any further 
discussion of them until Sect. IVII 
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FIG. 7: Graphs that contribute to a factorization anomaly. Other graphs that should be included are all other graphs with 
both gluons attaching the lower spectator to the and &4 lines at the upper part of the graph. In addition, the anomalous 
contribution to the TMD PDF of H2 is obtained by including the graphs with both gluons radiated from P2 and attaching the 
upper spectator to the k\ and k% lines. There is a total of 32 graphs of this type, including graphs with both gluons radiated 
from P2 . 



A. Graphs Contributing to Standard Factorization 

By repeating the steps of Sect. Mil graphs of the type 
shown in Figs. H] lead to the Wilson lines in Eqs. (fH)ll2U)) 
for the TMD PDFs of the standard TMD-factorization 
formula. The eikonal factors that arise from Figs. |U[a,b) 
correspond to keeping the order g Wilson line contri- 
bution from each of the TMD PDFs individually. Fig- 
ures |6]Jc,d) (and related graphs) correspond to keeping 
the order g 2 Wilson line contributions in one of the TMD 
PDFs, and the zeroth order from the other. 

The graphs considered up to now exhaust the contri- 
butions obtained by expanding each of the Wilson lines 
in the standard TMD-factorization formula Eq. (|24p up 
to order g 2 (in the Wilson line). Therefore, any re- 
maining uncanceled contributions from the graphs of the 
type shown in Figs. [JJ [8] violate the standard TMD- 
factorization formula. 



B. Violation of Standard TMD-Factorization 

The graphs of the type shown in Fig. where two 
gluons attach the spectator of one hadron to the op- 
posite side of the hard subgraph, have been shown in 
Refs. [ill, [2(| to be non- vanishing in calculations of both 
the unpolarized cross section and an SSA for an Abelian 
gauge theory. These graphs, therefore, yield a violation 
of the standard TMD-factorization formula Eq. (|24l) . 

The contribution to an anomalous SSA requires an 
imaginary part from the extra eikonal factors. In Ref. (T§j 
the imaginary part came from graphs with one gluon at- 
taching near to the parent hadron and the other attach- 
ing near to the opposite hadron (e.g., Fig. [^b)). These 



graphs vanish in our model because they again include a 
color factor Tr c [t a ] = 0. 

The non-cancellation of the unpolarized factorization 
anomaly, however, remains. It arises from graphs with 
both extra gluons radiating from the same hadron and 
attaching at the side of the hard subgraph near to the 
opposite hadron, as in Fig. [7J To make the violation of 
standard TMD-factorization explicit, we will briefly re- 
view the steps of Ref. [l8[ . We are for the moment only 
concerned with the unpolarized cross section so we will 
temporarily simplify the model of Sect. [II] by making all 
the hadron and spectator fields scalars. We further sim- 
plify the calculation by assuming m qi = m q2 — m^ 11 — 
m.0 2 = m q . 

For the graphs with the extra gluons on opposite sides 
of the cut (Fig.[7JJa) and related graphs), the factorization 
anomaly for the unpolarized cross section comes from 
replacing the extra off-shell quark propagators by their 
eikonal propagators. Including all ways of attaching £1 
and I2 to the &2 and lines, one finds the same result 
as in Ref. [l8j|, but with the Abelian <?i and gi charges 
replaced by the appropriate color factor: 

Tr c [t a t b ] <?X< + -^—) x 

\—li+ie l^ + iej 

\ — 12 — it I2 — ie) 

47ry«Tr c [t a t b ] S(l+)S(l+). (25) 

For the graphs with both extra gluons coupling on the 
same side of the cut (Fig. [TJb) and related graphs)), the 
steps are similar. Again including all ways of attaching 
l\ and I2 at the upper part of the graph, one finds 
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Tr c [t a t b ] 5 2 « 



1 \ / 1 



1 \ / 1 \ ( 1 \ / 1 



Tr c [i°t 6 ] 5 2 « f-^- + - J-^) f 7+^- + = - WnfnJ'&c [i * 6 ] 5{lt)5{l+). (26) 



The anomalous eikonal factor here is the same as in Eq. (|25j) . apart from an overall minus-sign. However, the 
propagator denominators for graphs like Fig. Eta) (both gluons on opposite sides of the cut) are different from the 
propagator denominators for graphs like Fig. [T^b) (gluons on the same side). For graphs with the extra gluons on 
opposite sides of the cut, one finds the following contribution to the TMD PDF of hadron H\: 

h{k 1T ) = 

g 2 A 2 Tr c [tH b ] Tr c [t b t a ] + f [2(p+ - fc+) + 1+] [2(p+ - fc+) + /+] 

(2tt)» XlPl J h (l\ + ^ e )(^ 2 - ie) [(ki hf - m§ + ie] [(h l 2 f -m\- ie] 

(2nrS(l+)5(l+)S((p 1 -k i r-mt) 



[(pi — fci + ^i) 2 - to 2 , + ie] [(p — fc + Z 2 ) 2 — "J 2 — «e 
ff 2 A 2 T 2 (7V 2 - 1)^(1 -X!) 



,. - l)3:i(l-3:i) f 2 2 -pr 1 

— — — d hra l 2T II , 2 r , — ro- 5T- ' 27 J 

2567T 7 7 j—12 l 0T k lT - ljT 2 + ™^ 



This is the same result as in Ref. [la], except that the gluon is massless and there is a non-Abelian color factor 
multiplying the integral. Equation (|26[) allows for a similar calculation of the remaining contribution to the TMD 
PDF for hadron H\ from the graphs with the extra gluons on the same side of the cut: 

T(k s -9 2 XjTUN!-l)x 1 (l-x 1 ) f 2 2, 1 



256^7 J q T q T [ (klT _ i 1T _ i 2T )2 + m 2] [ fc 



The mismatch in denominators between Eq. (|27j) and 
Eq. (|28|) means that the full contribution Ii(kix) + 
I%(kiT) does not generally vanish point-by-point in k\T- 
The sum of graphs like Fig. [7] therefore results in un- 
canceled terms that are not accounted for by the standard 
Wilson lines. Hence, the standard TMD-factorization 
formula Eq. ([24]) fails for unpolarized scattering. 

Exactly analogous observations apply to the TMD 
PDF of the other hadron if the two extra gluons are 
radiated from the spectator in H 2 , collinear to the mi- 
nus direction. In that case, the eikonal factors analo- 
gous to Eqs. (|25[ |2"6")) will instead use a vector n 2 and 
the delta functions from the eikonal factors will be 6(1^ ) 
and 5(1%) [33]. So, both TMD PDFs yield factorization 
anomalies at the two-gluon level in the unpolarized cross 
section. 



sequence of eikonal factors in the factorization anomaly 
terms of Eqs. (|25I26|) . including the trace around a color 
loop, require a color-traced Wilson loop operator to be 
inserted into the definition, Eq. (0, of the TMD PDF 
for Hi. Each of the eikonal factors in Eqs. (|25I26I) cor- 
responds to an attachment to a leg of the Wilson loop. 
Therefore, the standard TMD PDF in Eq. (J23]) for Hi 
should be replaced with, 

i[ni,(D)]/ , \ _ 
9 Hi \ x li k lT) — 



x ^ J 

x (Hi, si\4 r (0, w-,w t ) U% [0, w]U^ } 0i. s (O)|ffi, si). 

(29) 



V. GENERALIZED TMD-FACTORIZATION 



From the factorization anomaly terms Eqs. (1251126)) one 
may determine what the modified gauge link structure 
must be for the TMD PDF of hadron Hi in a general- 
ized factorization formula [9II12I [lEj. In our example, the 



This is the same as the standard TMD PDF definition in 
Eq. (|19|) apart from the insertion of the following color- 
traced Wilson loop operator: 

U^ = U^%w] (U n ^[0,w})ji = 
Tr c \v Q (ni)I(ni)V^(ni)V w (ni)I^(ni)V^(ni)\ . (30) 
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The sequence of eikonal factors in Eqs. (|25I26[) corre- 
spond exactly to gluons attaching to the "V" Wilson lines 
in Eq. (|30|) . As before, the transverse links at infinity do 
not contribute to leading power in Feynman gauge. The 
color trace corresponds to the trace over the color loop 
in the upper part of the graphs in Fig. [7] and is what 
is needed to get the factor of Tr c [t a t b ] in Eqs. ([251260 , 
Hence, the terms that violate standard factorization arise 
from expanding the operator in Eq. (|29[) up to order 

g 2 (see also Ref. 0). 

An exactly analogous TMD PDF is obtained for 
hadron H 2 when both extra gluons are radiated from H 2 
in the collinear minus direction and attach to the lower 
half of the graph. The main legs of the Wilson loop in 
the TMD PDF for H 2 point in the direction of a vector 
n 2 : 



< 2 ' (D) Wfc2T) 



X2P2 



dw + d 2 xv t 
(2tt) 3 ' 



-iX2Pi iu + +*kfw f 



x (H 2 , s 2 \4 r (0, w + , w t ) C/; i s 2 [0, w]Ufa 4>2,M\H2, s 2 >. 

(31) 



Thus, the factorization anomaly terms from Sect. II V 51 
specify the type of new Wilson line insertions that are 
needed if generalized TMD-factorization is to hold. The 
generalized TMD-factorization formula is: 

da In ® $ [ £>™(x u k 1T ) ® ^l' [a) \x 2l k 2T )® 

® 5 [2) (ki T + k 2T - k 3 T - k 4T ). (32) 

Again, % is the same hard factor that appeared at zeroth 
order in Sect. [HI The only difference from Eq. (|24|) comes 
from the Ujj-^ insertions in the TMD PDFs, indicated by 
the (□) superscripts on the individual TMD PDFs in 
Eq. (|3"2"1) . We have again included a "!" over the equal 
sign to indicate that even this formula will ultimately fail 
when the remaining graphs are considered. 

Equation fl32J|, with the TMD PDFs defined as in 
Eq. (|2T))) , follows the method proposed in [9|, [12J, ll5j , and 
is the most natural generalization of TMD-factorization. 
For the purpose of our counter-argument to generalized 
TMD-factorization, the crucial point is that the new Wil- 
son line insertion Eq. (fBTJf involves a trace over color. 
If we instead tried contracting the color indices of the 
new Wilson lines with the quark fields at points and 
£, then we would obtain the wrong color factor for the 
anomalous contribution in Eqs. (|25I26[) . An even more 
serious problem would be that the generalized TMD- 
factorization formula would produce spurious terms at 
the one-gluon level, already treated in Sect. IIIII Namely, 
single-gluon attachments to the extra legs of the Wilson 
line would introduce factorization anomaly contributions 
analogous to Eq. flTTl) . but with different non- vanishing 
color factors, contradicting the analysis from Sect. IIIII 
which found no violation of standard TMD-factorization 



at order g in the expansion of the Wilson lines (polarized 
or unpolarized). Thus, consistency between the order- 
g and order-<? 2 contributions to the Wilson lines requires 
that the new Wilson line structures include the trace over 
color as in Eq. (|30)) . 

Other possible redefinitions of the TMD PDFs are un- 
related to the color structure of the Wilson lines in the 
main correlation functions. These include the modifica- 
tions mentioned in Sect. HTJ such as tilting the direction 
of the Wilson line slightly away from the light-like direc- 
tion. Thus, accounting for all contributions up to order 
g 2 , in the Wilson line, and requiring consistency between 
different orders in <?, strongly constrains the possibilities 
for a generalization of TMD-factorization. 

The order-<? contribution from a Wilson loop in 
Eq. (|32|) is zero simply because, with one gluon attach- 
ing to the Wilson loop, the color trace in the definition of 
the Wilson line operator Eq. (|3Tfl) involves only one color 
generator Tr c [t a ] = 0. So, 



rrni 
U (D) 



= U n2 

O(g) (D) O(g) 



0. 



(33) 



Therefore, there is no disagreement between Eq. 
and standard factorization, Eq. (|24|) . from graphs with 
only one extra gluon. To find a disagreement between 
the standard and generalized TMD-factorization formu- 
las (Eqs. JM]) an d (|3"2")1 ), at least two gluons need to be 
collinear to one of the hadrons as in Sect. IIV Bl 



VI. BREAKDOWN OF GENERALIZED 
TMD-FACTORIZATION 

We will now directly illustrate a breakdown of the gen- 
eralized TMD-factorization formula Eq. ([32]) by calculat- 
ing the anomalous two-gluon contribution to a double 
Sivers effect from the graphs in Fig. [5] First, let us sum- 
marize the situation so far: 

• Considering graphs with up to two extra glu- 
ons, there is a one-to-one correspondence between 
graphs of the type shown in Fig. [3j/Fig. [6] and the 
contributions to the Wilson lines in the standard 
TMD-factorization formula Eq. (|2"i|). 

• The sum of graphs with two extra gluons radi- 
ated from one spectator and attaching on the oppo- 
site side of the hard subgraph (as in the graphs of 
Fig- EJ) results in terms that violate standard TMD- 
factorization in unpolarized scattering. These are 
the contributions already discussed in Ref. [l8j . 

• The two-gluon factorization anomaly terms found 
in Sect. IIV Bl specify which modifications of the 
Wilson lines in the TMD PDFs are needed if a form 
of factorization is to be recovered. It is found that 
the modified TMD PDFs must each contain an ex- 
tra color-traced Wilson loop, as in Eq. (j2"5)) . This 
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FIG. 8: Graphs that contribute to a violation of generalized TMD-factorization. Other graphs that should be included are 
those with all possible attachments of l\ to the k± and lines, and all possible attachments of h to the k^ and ki lines, and 
all Hermitian conjugate graphs. In total there are 16 graphs of this type. 



result also follows from the same steps for find- 
ing general Wilson line structures from low order 
graphs as discussed in Rcf. 12]. 

• Because of Eq. (p£?)) . the generalized TMD- 
factorizaton formula Eq. (|32|) cannot include any 
other two-gluon contributions to Wilson lines. In 
particular, there can be no contribution that cor- 
responds to graphs like Fig. [51 where one gluon is 
radiated collinear to each hadron simultaneously. 
These graphs would correspond to separate order- 
s' contributions from the Wilson lines in separate 
TMD PDFs. If they are non-zero, then they 
contribute to a violation of both standard TMD- 
factorization in Eq. (f2~4"|) and the generalized TMD 
factorization formula in Eq. (|32[) . To incorporate 
such a contribution, one would have to modify the 
Wilson line in each TMD PDF such that it in- 
cludes a single-gluon contribution to a factorization 
anomaly. But this would contradict Sect. Mil where 
it was shown that there is no violation of standard 
factorization with just one gluon. Hence, contribu- 
tions from graphs like Fig. [5] cannot be consistently 
incorporated into a generalization of factorization 
simply by modifying Wilson lines in separate corre- 
lation functions. If they give a non-vanishing con- 
tribution, then there is a clear violation of general- 
ized TMD-factorization. 

We will therefore prove that generalized TMD- 
factorization, Eq. (j3"2")l . is violated by showing that the 
sum of graphs of the type illustrated in Fig. [5] give a 
non- vanishing contribution to a DSA. 

First, we note that all graphs of the type shown in 
Fig. [8J include the non-zero color factor 



Tr c [t a t h ] Tr c [t b t a ] 



T 2 {N 2 C - 1). 



(34) 



Next, we must ensure that there is no cancellation be- 
tween graphs. 



A. Same Side of the Cut 

In the sum of graphs like Fig. Eta), where both gluons 
are on the same side of the cut, the eikonal factors give 
a total contribution equal to 



1 



1 



1 



1 



—li+iej \—l 2 + if J Vi + if J \— Z 2 + if 
1 \ / 1 \ ( 1 \ ( 1 



l~i + it J \l 2 + if J \l\ + if J \Z 2 + if 

= -*n*8(l+)8(i;). (35) 

Since spin dependence is needed in both H\ and H 2 for 
a DSA, then there are also two factors corresponding to 
Eq. (|9]) (but with one corresponding to a p 2 spectator 
attachment). Taking into account both of the resulting 
factors of i gives an overall factor of i 2 = —1. Combined 
with Eq. (|35[) . the relevant factor from extra collinear 
gluons is then 4tt 2 5(1+)8(1 2 ). The same result is obtained 
from the Hermitian conjugate graphs. 



B. Opposite Side of Cut 

The sum of graphs with one extra gluon on each side 
of the cut, as in Fig. [8jb), works in much the same way. 
The eikonal factors give 



1 



1 



-Zj — if ) \— 1 2 + if J \Zj — if J \— 1 2 + if 



1 



1 



— Zj — ie J \l 2 + if J \li — if J \l 2 + if 

^4tt 2 6(1+)5(1 2 ). (36) 



For a DSA, there is a factor of i from a factor analogous 
to Eq. © for the pi-spectator attachment of the gluon 
on the left side of the cut and a factor of — i from a 
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factor analogous to Eq. (fTTj) (but for a gluon attaching 
at the ^-spectator) on the right side of the cut, giving an 
overall factor of i{— i) — +1. So, combined with Eq. (|36|) 
the relevant factor from extra gluon attachments is again 



itly, one can use Eqs. (l36|) (extracting the overall fac- 
tor of i(— i) = 1 that comes with the two spectator at- 
tachments) to explicitly calculate the contribution from 
graphs with gluons on opposite sides of the cut: 



C. Together 

Summing all graphs of the type shown in Fig. [51 there- 
fore, results in just a single integral. To check it explic- 



2T 2 {N 2 - 1)<? 4 f d 4 h f d 4 k 2 f d% 



2s J (2tt) 4 J (2tt) 4 J (2tt) 4 J (2tt) 



vi 



(27r) 4 <5 4 (fci + k 2 - h - k 4 )(2Tr) 2 5(l+)S(l 2 -)x 



,2,,2(fci + k 3 ) ■ (k 2 + k 4 )\ 2 2ej k s\l\-p\ (m Hl { 1 ~ x i) 



Al ^ (fci - fc 3 ) 2 - M\ J [l\ + ie][(h -hf- ml + ie][(pi - h + h) 2 - m 2 h + ie][k 2 

2ey k 's{l 2 p 2 {m Ho (l - x 2 ) + m^ 2 ) 2 2 2 2 2 

(37) 

Here we have dropped the irrelevant factors of (2E 3 ^), Xi(2), and (2ir) 6 which appeared in Eq. We have also 
dropped the l\ and l 2 dependence inside the factor in braces. This is permitted because the 5-functions will set if 
and l 2 to zero, and because the other components yield power suppressed corrections. The 8- functions can be used 
to evaluate the kf^k^jk^, k^ , if, and k 2 y integrals, and the and l 2 integrals can be evaluated by contour 
integration. The result is 



T 2 (N 2 - 1) 5 4 



■H(fci,fc 3 ,fc 4 )x 



(2tt) 2 J (2tt) 2 J (2tt) 2 

[tjks{l\{m Hl {\ - xi) + m^_)[[ej^s| 1%' (m H2 (l - x 2 ) + m fe )| 

iir f 2r[*iT + • • • ] [Oir - k 1T ) 2 + ■ ■ • ] [(k 3T + k 4T - k 1T ) 2 + ■ • • ] [(1 2T - k 3T - k 4T + k ir ) 2 + 



(38) 



The factor H(ki, fc 3 , ki) is now the hard part in Eq. ([6]). 
The symbol "• ■ ■ " in the denominators refers to terms 
that involve only x%, x 2 and masses. For graphs with 
both gluons on the same side of the cut, exactly the same 
result (with the same overall sign) is obtained by using 
Eq. (|35|) and taking into account the overall factor of 
i 2 = — 1 that comes from the two spectator attachments. 
There are, again, identical results coming from the Her- 
mitian conjugate graphs. 

The integral in Eq. ([38]) is not generally zero as can be 
checked by considering a hxed q = k 3 ^ + k 4 x. How- 
ever, there is no corresponding factorization anomaly 
term contained in an expansion of Eq. (j32j) up to or- 
der g 2 in the Wilson lines. The modified Wilson lines 
in the TMD PDFs of the generalized TMD-factorization 
formula are specified by the factorization anomaly terms 
that arise when gluons are radiated from just one hadron 
at a time, as in Sect. II VI However, the resulting general- 
ized TMD-factorization formula is in contradiction with 



graphs like Fig. [5J where there are simultaneously glu- 
ons radiated from Hi and H 2 . Thus, generalized TMD- 
factorization breaks down at the two-gluon level. 

Indeed, at up to two gluons, the generalized TMD- 
factorization formula Eq. Q32p gives no contribution to a 
DSA beyond what is predicted by the standard TMD- 
factorization formula Eq. (|24p . The only factorization 
anomaly contribution from Eq. (|32l) is to the unpolar- 
ized cross section, and arises from graphs with both extra 
gluons radiated from the same hadron. There must be 
at least one gluon coming from each hadron at the same 
time to have spin dependence in both hadrons simulta- 
neously. So, the only contribution to a double Sivers 
effect that violates standard TMD-factorization is from 
the sum of graphs in Fig. [5J and these are not accounted 
for in Eq. (|32|) . To recover the factorization anomaly in 
the DSA, one would have to somehow allow each TMD 
PDF to have a non-vanishing factorization anomaly con- 
tribution at the level of one extra gluon. This, however, 



16 



would contradict the single-extra-gluon treatment of the 
SSA in Sect. IIIII Therefore, a TMD-factorization formula 
cannot be recovered by simply modifying the Wilson lines 
in each TMD PDF separately. We have thus illustrated 
a specific example of a violation of generalized TMD- 
factorization. 

One can interpret the breakdown of generalized TMD- 
factorization visually by imagining the flow of color in 
the factorized expression, Eq. as compared to the 

original unfactorized graph. The graph in Fig. [8J5, for 
example, must correspond to the product of two single- 
gluon contributions to the Wilson lines in Eq. (1521) , if that 
formula is correct. The resulting color structure can be 
visualized as in Fig. [9] The narrow double lines represent 
the Wilson lines in each TMD PDF. The box shaped Wil- 
son lines represent the extra color-traced Wilson loops in 
Eq. The color trace over each of the Wilson loops 

means that only a color singlet gluon can be exchanged 
between a Wilson loop and a spectator at the single-gluon 
level, as illustrated by the factors representing the TMD 
PDFs in Fig. [9] Therefore, these contributions to the 
TMD PDFs must vanish. By contrast, in the original 
unfactorized graph, the interlaced color flow of the two 
gluons means there is no problem with exchanging non- 
singlet color, as illustrated by the color flow diagram in 
Fig. [TU] This corresponds to the non-zero color factor in 
Eq. 



VII. DISCUSSION AND CONCLUSION 

We have shown that a generalized TMD-factorization 
formula like Eq. Q /Eq. (|32p does not in general exist for 
back-to-back hadro-production of high-p t hadrons due 
to the non-trivial interplay of gluons that are collinear 
to both hadron directions. The failure of generalized 
TMD-factorization calls into question results that use 
it as a starting assumption. This includes calculations 
of weighted cross sections which retain a memory of the 
non-universality of the TMD PDFs. 

It is interesting to note that the same complications 
do not seem to arise if the gauge field is Abelian. In 
that case, the Abelian "color factors" can be identi- 
fied with coupling constants, such as in the model from 
Ref. [17j . Then the contribution to a DSA from graphs 
like Fig. |8] is indeed just a product of single gluon contri- 
butions to Wilson loops from separate correlation func- 
tions (which are non-zero in the Abelian case). There- 
fore, there is no obvious contradiction with generalized 
TMD-factorization for the Abelian case, at least at the 
level of graphs discussed here. While a full proof of gen- 
eralized TMD-factorization for the Abelian gauge theory 
does not yet exist, the breakdown of generalized TMD- 
factorization that we have found is a specific consequence 
of the non- Abelian nature of QCD. 

It is worth stressing that the failure of generalized 
TMD-factorization occurs in a regime where factoriza- 
tion would ordinarily be expected to apply; namely for 



hard processes with large-p*. Moreover, it should not 
be thought that the breakdown of generalized TMD- 
factorization is specific to a DSA or SSA. We have calcu- 
lated for a DSA with color singlet boson exchange for the 
hard part in order to maximize the simplicity of the ar- 
gument, but the reasons for the breakdown are quite gen- 
eral. The counter-example we have provided in this paper 
is sufficient to show that a general proof does not exist. 
That is, it is not possible to consistently define TMD 
PDFs for each hadron separately, even if we allow for 
process dependence in the Wilson line structures. To en- 
counter the same complications with generalized TMD- 
factorization in an unpolarized cross section, one needs 
the type of non-cancellation that was found in Sect. lIVBl 
(which comes from having two gluons collinear to one 
hadron) , but with at least one more gluon collinear to the 
opposite hadron. An example of the type of graph which 
can produce a violation of generalized TMD-factorization 
in the unpolarized cross section is shown in Fig. 1111 

Furthermore, the contribution from extra gluon at- 
tachments (which lead to factorization breaking) should 
not be thought of as negligible higher order corrections. 
They correspond to soft and collinear divergences in 
higher order hard scattering calculations. In real QCD, 
they are non-perturbative gluons in the strong coupling 
regime. For the perturbation series to have sensible con- 
vergence properties, it must be possible to rearrange 
terms such that the extra soft and collinear gluons are 
resummed to all orders into TMD PDFs. We also remark 
that, at the level of many extra gluons, there are more 
graphs with gluons collinear to both hadrons simultane- 
ously (which led to a breakdown of generalized TMD- 
factorization) than graphs with gluons radiated from one 
hadron only. 

While the observation of generalized TMD factoriza- 
tion breaking leads to frustrating practical difficulties 
in cross section calculations, it should not necessarily 
be regarded as a purely negative result. The question 
of whether transverse momentum effects can be mean- 
ingfully associated with parton transverse momentum in 
separate parton correlation functions for each hadron is 
intrinsically important in the search for an improved fun- 
damental understanding of QCD dynamics in hard col- 
lisions. A counter-proof implies the existence of effects 
which challenge normal partonic intuition, and suggests 
new avenues of research. While naive factorization fails, 
the fact that the extra gluons eikonalize suggests that 
Wilson lines still may play a role. Insight might be 
gained, for example, from methods currently being ap- 
plied to small- x physics (e.g. [H, HI])- Another possi- 
bility is to model factorization breaking effects by di- 
rectly calculating factorization breaking phase contribu- 
tions in perturbation theory, but with explicit infrared 
and collinear cutoffs. 

It is possible to understand the origin of the general- 
ized TMD-factorization breakdown intuitively as arising 
from non-linear effects in the phases acquired by partons 
as they pass through the A + and A~ fields of the collid- 
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da ~ W 





FIG. 9: Color flow resulting from the single gluon contributions for each of the Wilson loops in the TMD-factorization formula 
Eq. (|32p . H is the standard zeroth order hard part and the second two factors are the TMD PDFs. The narrow double lines 
represent Wilson lines. The boxes associated with each of the TMD PDFs correspond to the Wilson loops. The thick solid red 
and dotted blue lines (color online) illustrate the flow of color in each TMD PDF. Each of the contributions to a TMD PDF 
shown here is exactly zero because each includes a factor Trc [t a ] = 0. 
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FIG. 10: Color flow in the unfactorized graph with a single gluon collinear to each of the incoming hadrons as in Fig. [8] The 
thick solid red and dotted blue lines (color online) again illustrate the flow of color. Non-singlet color can easily be exchanged 
and results in a non-zero contribution. Compare with Fig. [9] 



ing hadrons. If the overall phase were simply the product 
of the phases induced by the A + fields from hadron Hi 
and the A~ fields from hadron H 2 , then one could asso- 
ciate any process-dependent phases induced by the A + 
field in hadron Hi with a modified Wilson line for the 
TMD PDF of Hi and, likewise, any process-dependent 
phases induced by the A~ field from H 2 could be asso- 
ciated with a modified Wilson line for the TMD PDF of 
H 2 . However, in the non-Abelian theory the role of the 
A" gluons in H 2 is affected by the presence of the A + 
gluons from Hi and visa- versa. A direct example of this 
is Fig. |8j/Eq. (l38l) . where a single A~ gluon exchanged 
between H 2 and the opposite-side struck quark gives a 
non-zero contribution, but only because there is simul- 
taneously an A + gluon exchanged between Hi and the 
other struck quark. This means that one cannot address 
the role of phases induced by the A + and A~ fields inde- 



pendently, but instead must deal with them simultane- 
ously. The result is a kind of nonperturbative correlation 
which cannot be identified as arising strictly from gluons 
coming from either hadron independently, but only from 
the combination. 
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Because in Eq. Q there is an integration over all /cit, 
one generally must be careful in treating the region where 
the t goes to zero in the hard part. For this paper, we are 
specifically interested in the collinear region where kir 
is small, so we will assume that any contribution from t 
close to zero is removed by a cutoff on large kir- 
In a general proof of factorization, one may need to worry 
about cancellations in the "supersoft" region, kir « 
Aqcd, since we are now using a massless theory, in con- 
trast to Ref. [l7)| . In this paper, we are concerned mainly 
with gluons with kir ~ Aqcd so these issues do not 
affect our argument. 

In earlier sections we neglected gluon exchanges between 
active quarks because they do not lead to spin depen- 
dence. In the unpolarized cross section, the complete re- 
sult is obtained only after including interactions between 
active quarks. In the collinear regions they result in the 
same set of eikonal factors as in Eqs. (1251 1261) and anal- 
ogous Wilson line contributions. The key point is that 
the usual steps that normally would lead to a cancella- 
tion of the factorization anomaly between different cuts 
of the same graph in collinear factorization, fails in the 
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transverse momentum dependent case. 



